Systematic behaviour of the in-plane penetration depth in d-wave cuprates 
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We report the temperature T and oxygen concentration dependences of the penetration depth of 
grain-aligned YBa2Cua07_i with S= 0.0, 0.3 and 0.43. The values of the in-plane \ a b{0) and out- 
of-plane A c (0) penetration depths, the low temperature linear term in A a i>(T), and the ratio A c (0) 
/Aab(T) were found to increase with increasing 5. The systematic changes of the linear term in 
A a (,(T) with T c found here and in recent work on HgBa2Ca n _iCu n 02 n +2+,5 (n = 1 and 3) 1 are 
discussed. 



In a recent study 1 of the c-axis coupling of d-wave 
high-T c cuprates we reported the values and tempera- 
ture (T) dependences of the in-plane (X a b) and out-of- 
plane (A c ) penetration depths for slightly overdoped 2,3 
HgBa2Cu04 + 5 (Hg-1201) with critical temperature T c 
— 93 K and slightly underdoped 4 ' 5 HgBa 2 Ca2Cu308+5 
(Hg-1223) with T c = 135 K. For both compounds the 
low temperature dependence of A a & was found to be lin- 
ear as expected for d-wave superconductivity. In fact 
normalised plots of [A a b(0) / A a h(T)] 2 versus T/T c were 
the same, and like YBa2Cu307 (YBCO7) 6 agreed very 
well with mean field (MF) theory for a weak coupling 
d-wave superconductor. However, recent angle-resolved 
photoemission spectroscopy (ARPES) 7 and tunneling 8 
data strongly suggest that for underdoped samples the 
superconducting gap, (Ao), remains constant, or even in- 
creases slightly, while T c falls and so large deviations from 
MF theory might be expected. We have therefore ex- 
tended our investigation to deoxygenated (underdoped) 
pure YBa2Cu307_a using the same ac susceptibility tech- 
nique to measure the penetration depth. 1 ' 6,9 ' 10 

We report experimental results for the values and tem- 
perature dependences of A Q & and A c of high quality c - axis 
grain-aligned orthorhombic 11 YBa2Cu307_^ (which has 
two CuC>2 planes per unit cell as well as Cu-0 chains) 
with 5 = 0.0, 0.3 and 0.43, and compare them with 
tetragonal 12 Hg-1201 with one Cu02 plane per unit cell 
and tetragonal 13 Hg-1223 with three CuC>2 planes per 
unit cell. We find that the presence of the linear term 
in A Q &(T) is independent of the number of Cu02 planes 
per unit cell, carrier concentration, crystal structure, 
anisotropy and the presence of chains. Surprisingly our 
data show good agreement with weak coupling d- wave 
theory, and the linear term in [A a &(T) / A a &(0)] appears 
to scale with T jT c . This result highlights the need for 
detailed consideration of the relationship between super- 
conducting and normal state energy gaps in underdoped 
cuprates. 

Sample preparation was carried out by the standard 
solid state reaction process using high purity (99.999%) 
Y2O3, BaCC-3 and CuO oxides. Electron probe micro- 
analysis and x-ray diffraction showed that all samples 
were single phase within an accuracy of rs 1%. The fully 
oxygenated, S = 0.0 (YBCO7, T c = 92 K), samples were 
prepared by annealing bulk pieces in pure oxygen atmo- 



sphere at 380 °C for 24 hours and then slowly cooling to 
room temperature. (Hereafter T c represents the temper- 
ature where the onset of superconductivity occurs in the 
ac susceptibility data for a measuring field H ac = 3 G 
rms and frequency f = 333 Hz.) The 5 = 0.3 (YBC0 6 . 7 , 
T c = 66 K) sample was prepared by annealing in pure 
oxygen atmosphere at 650 °C for 12 h and then quench- 
ing in liquid nitrogen, while the 5 — 0.43 (YBCO6.57, T c 
— 56 K) sample was annealed in 0.2% O2 / N2 atmo- 
sphere at 550 °C for 12 hours and also quenched into liq- 
uid nitrogen. The final oxygen contents were determined 
from the weight change of a fully oxygenated reference 
sample. The 5 = 0.0 bulk piece was lightly ground and 
sedimented in acetone to obtain a well-defined grain size 
distribution. The sedimented powders were then heat 
treated to repair any structural damages to the surface 
of the grains. 14 For 5 = 0.3 and 0.43 on the other hand 
a bulk piece for each d was lightly ground and sieved 
through a 20 (an sieve in an argon glove box to obtain a 
well-defined grain size distribution 10 and to avoid surface 
degradation of the crystallites. 14 The collected powders 
were then kept in argon atmosphere for 30 min before be- 
ing aligned. All powders, 5 = 0.0, 0.3 and 0.43, were mag- 
netically aligned in epoxy as described earlier. 1 ' 6,10 The 
average grain diameters corresponding to the 50% cumu- 
lative volume point were 5 and 10 /jm for the fully oxy- 
genated and the oxygen deficient samples, respectively. 
The fraction of the unoriented powder in all grain aligned 
samples was estimated to be < 5 %. Rocking curve anal- 
ysis of the S = 0.0 and S > 0.0 samples gave a full width 
at half maximum of ± 1.4 and 1 °, respectively. 15 Low 
field susceptibility, Xj measurements were performed us- 
ing commercial equipment (down to 4.2 K) for samples 
with S = 0.0, 0.3 and 0.43. The sample with S = 0.43 was 
also measured down to 1.2 K using a home built suscep- 
tometer. Details of the experimental technique and the 
application of London's model for deriving A from the 
measured x i n cuprate superconductors can be found in 
earlier publications. 1,6 ' 9,10,16 

The values of A a f,(0) derived from our data are 0.14, 
0.21 and 0.28 /im and the corresponding values for A c (0) 
are 1.26, 4.53 and 7.17 ^m for S = 0.0, 0.3 and 0.43, 
respectively. The errors in A a h(0) arising from a pos- 
sible uncertainty of ± 5% in the alignment can be as 
high as ± 25%, whereas those in A c (0) are ± 8%. How- 



ever, the corresponding uncertainty in the linear term in 
[A o b(T)/A a fc(0)] is much less, at most ± 10%. The present 
results differ from previous work 17 in which the surfaces 
of the particles were probably not as clean and the degree 
of grain alignment was probably lower. As T c is reduced 
by lowering the carrier concentration (for S = 0.3 and 
0.43), [1 / A a b 2 (0)] falls, a behaviour which has been ex- 
tensively discussed in terms of the Uemura relation. 18,19 
The ratio 7 = [A c (0) / A a b(0)] i.e. the anisotropy, in- 
creases with oxygen deficiency. 

Figures 1(a) and 1(b), show characteristic low temper- 
ature plots of [A(T) / A(0)] for the ab- plane (measured 
with the applied field H || c) and c-axis (measured with 
H || ab), respectively, for the three oxygen concentrations 
studied. The low temperature (T / T c < 0.25) linear term 
in A a t(T), is 4.8 A/K for YBCO7 in good agreement with 
that found from microwave measurements on YBCO6.95 
single crystals. 20 As oxygen is removed from the lattice 
(the chains) the linear term increases to 12 and 20 A /K 
for 5 = 0.3 and 0.43, respectively. For YBCO7 we also 
observe a linear T dependence in A c at low temperatures 
but the relative change is about a factor of two smaller 
than in [A a6 (T) / A a6 (0)], while A C (T) of YBC0 6 . 7 and 
YBCO6.57 obeys a T 2 behaviour at low T. Details of the 
systematics of A C (T) of cuprate superconductors can be 
found in Refs [1,6,21,22]. 

In Fig. 2 we present normalised plots of [A a {,(0) 
/A a6 (T)] 2 and [A c (0) /A C (T)] 2 versus T /T c . There is ex- 
cellent agreement between the [A a fc(0) / A a t(T)] 2 curves 
for the three oxygen concentrations. The data in Fig. 
2(a) are compared with the weak coupling theory for a 
d-wave superconductor (solid line). 23 It can be seen that 
the d-wave curve fits the data very well. On the other 
hand the [A c (0) / A C (T)] 2 curves do not fit the d-wave 
curve and also differ from each other slightly, because of 
the effect of the interplane coupling on Ac(T). 1 ' 6 ' 10 ' 21 ' 22 
We also find that [A afc (0) / A a6 (T)] 2 of YBa 2 Cu 3 7 -«5 
agrees with that of Hg-1201 and Hg-1223. 1 ' 10 The be- 
haviour of [A a {,(0) / A a h(T)] 2 is generally similar to that 
of YBa 2 (Cui_ x Zn x ) 3 7 (x =0.02 and 0.03), 6 except at 
very low temperatures where a T 2 term developed in the 
Zn doped samples due to impurity scattering. 

The full temperature dependences of [A a &(0) / 
A a6 (T)] 2 for YBa 2 Cu 3 7 -«5 (5 = 0.0, 0.3 and 0.43), 
YBa 2 (Cui_ x Zn x ) 3 7 (x =0.02 and 0.03), Hg-1201 and 
Hg-1223 are also in agreement with recent data of 
Bi 2 Sr 2 CaCu 2 8+(5 (Bi-2212) 24 and Tl 2 Ba 2 Cu0 6+(5 (Tl- 
220 1) 25 single crystals measured by a microwave tech- 
nique with H || c, but they only agree with another set of 
microwave data (H || c) for Bi-2212 single crystals 26 at T 
/T c < 0.3. At higher temperatures the data in Ref. [26] 
deviate from the weak-coupling d-wave calculation. In- 
dependent evidence for the scaling behaviour of [A a & (0) / 
A a {,(T)] 2 with T /T c can also be found in a recent publica- 
tion by Bonn et al, 27 who measured the relative changes 
in A with temperature for underdoped, optimally doped 
and slightly overdoped untwinned YBa 2 Cu 3 7 _,5 crystals 
using a microwave technique and H _L c. However, in Ref. 
[27] the changes of [A„, 6 (0) / A a , fc (T)] 2 for YBC0 6 . 95 , at 



high temperatures, and YBCO6.6, over the whole tem- 
perature range, were smaller than ours and closer, at 
high temperatures, to the Bi-2212 data in Ref. [26]. We 
do not know the precise origin of this difference but we 
believe that for weakly coupled layers, data taken with 
H || c give the best measure of the superfluid density. 

Figure 3 shows plots of l-[A afc (0) / A a6 (T)] 2 vs T which 
is equivalent to [n s (0) - n s (T)] / n s (0), i.e. the nor- 
malised density of quasiparticle excitations, where n,,(T) 
is the density of condensed electrons at a temperature 
T, for all the samples studied. It is clear that the linear 
term in l-[A a b(0) / A a ;,(T)] 2 increases as T c decreases. 
If we use the standard BCS result for A(T) of a d-wave 
superconductor, 28 
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to fit the experimental data shown in Fig. 3 (at T / 
T c < 0.25) we find that Ao scales approximately with T c 
[Fig. 3 (inset)], giving A ~ 2T C , a value close to that ex- 
pected for weak- coupling superconductivity. 28 For com- 
parison we also include data for the s-wave perovskite 
Ba . 6 K .4BiO 3 (BKBO). 10 The compounds Bi-2212 24 ' 26 
and Tl-2201 25 would also give A « 2T C on this plot. 
The maximum error in the linear terms i.e. the values of 
A in Fig. 3 (inset), is ± 20 %. 

The scaling of Ao with T c , Fig. 3 (inset), is in 
agreement with early tunneling spectroscopy data 29 for 
several cuprates as a function of carrier concentration, 
ranging from the underdoped to the optimally doped 
regimes. It is not consistent however, with more recent 
tunneling 8 and ARPES 7 results for underdoped cuprates 
where Ao was actually found to increase slightly while 
T c falls. There seems to be two possible ways of ac- 
counting this discrepancy. One is that the recent spectro- 
scopic experiments 7,8 actually measure the normal state 
gap. In this scenario the effect of the normal state gap 
would be to leave small pockets of holes whose supercon- 
ducting properties are still described reasonably well by 
MF theory. The other is similar to a recent phenomeno- 
logical approach. 30 As shown in Fig. 4, it is probable, 
that within experimental error, the unnormalised plots of 
[l/A a h(T) 2 ], i.e. n s (T), versus T at low temperatures are 
parallel for samples with different T c values. This would 
correspond to the same number of excited quasiparticlcs, 
[n s (0) - n s (T)], at a given temperature for all T c values- 
as implied by specific heat work on underdoped YBCO. 31 
Such parallel shifts give [n s (0) - n s (T)] = aT, where a is 
independent of doping level (T c ). In combination with 
the well-known Uemura relation n s (0) oc T c , 18 this gives 
[1- n s (T) /n s (0)] = b T /T c , where b is independent of 
T c . So at low T, [A ah (0) / A ah (T)] 2 versus T/T c would 
still scale on to a single curve even when the MF relation 
Ao / T c « 2, is strongly violated. 

In conclusion, we have studied A a t(T) and A C (T) of 
high quality grain-aligned YBa 2 Cu 3 7 _,5 with S = 0.0, 
0.3 and 0.43. The values of A a &(0), A c (0) and 7 were 
found to increase with oxygen deficiency. We find that 
the existence of the linear term in A a h(T) is independent 



of the number of CuC>2 planes per unit cell, carrier con- 
centration, crystal structure, anisotropy and the presence 
of chains. If viewed in isolation, all the penetration depth 
data presented here and most of the microwave measure- 
ments for H || c appear to be in excellent agreement with 
mean field theory for a weak coupling d-wave supercon- 
ductor for which Ao / T c w 2. However, recent spectro- 
scopic data are more consistent with a different approach 
30 in which there is a strong interplay between the super- 
conducting and normal state gaps. Clearly the relation- 
ship between these two gaps is of crucial importance for 
understanding superconductivity in the cuprates. 
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FIG. 1. Low temperature plots of (a) [A a f>(T) / A o f>(0)] 
and (b) [A C (T) / Ac(0)] for YBC0 7 (closed circles), YBC0 6 . 7 
(open circles) and YBCO6.57 (open squares). The T c , A o t(0) 
and A c (0) values are given in the text. 
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FIG. 2. Plots of (a) [A a6 (0) / \ ah (T)} 2 and (b) [A c (0) / 
A C (T)] 2 as functions of T /T c for YBC0 7 (closed circles), 
YBCOe.7 (open circles) and YBCOe.57 (open squares). The 
solid line in (a) is the theoretical prediction for the normalised 
superfluid density from the weak coupling BCS theory for a 
d-wave superconductor. 23 
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FIG. 4. Low temperature plot of [l/A a(j (T)] (i.e. n 3 ) for 
YBCO6.57 (open squares) and Hg- 1201 (open triangles) show- 
ing the approximate parallel shift of n s with T c as discussed 
in the text. The dashed lines, immediately above and be- 
low each data set, indicate the maximum possible error in 
1/A 2 6 (T) arising from ± 5% uncertainty in the alignment (see 
text). 
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FIG. 3. Low temperature plot of 1- [\ ab (Q) / A a6 (T)] 2 for 
Hg-1223 [A a6 (0) = 1770 ± 300 A] (closed triangles), 1 Hg-1201 
[Aab(0) 1710 ± 250 A] (open triangles), 1 YBCO7 (closed cir- 
cles), YBCO6.7 (open circles) and YBCO6.57 (open squares). 
Inset: T c versus 2Ao as derived from the plot in the main 
panel (see text for details). BKBO is included for compari- 
son. The solid line is drawn as a guide to the eye. 



